Complex Phase Behavior of The System of Particles with Smooth Potential with 

Repulsive Shoulder and Attractive Well 
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We report a detailed simulation study of the phase behavior of core softened system with attractive 
well. Different repulsive shoulder widthes and attractive well depthes are considered which allows 
to monitor the influence of repulsive and attractive forces on the phase diagram of the system. 
Thermodynamic anomalies in the systems are also studied. It is shown that the diffusion anomaly 
is stabilized by small attraction. 

PACS numbers: 61.20.Gy, 61.20.Ne, 64.60.Kw 



I. INTRODUCTION 

In 1970-th Hemmer and Stell introduced a model sys- 
tem which can effectively approximate real substances 
0, Q . The model they introduced belongs to the class of 
potentials which nowadays is called core-softened poten- 
tials, i.e. the system with softening of the repulsive part 
of interparticle interaction. Hemmer and Stell showed 
that such systems can demonstrate some unusual phase 
transitions - transition between two crystal phases of the 
same symmetry [l|, Q which ends in a critical point sim- 
ilar to the gas-liquid one, and liquid - liquid phase tran- 
sition. In Ref. [3| it was proposed that the high density 
transition with critical point qualitatively represents the 
isostructural transition in Ce and Cs. 

Great interest to this kind of systems appeared after 
Hemmer and Stell publications. Many researchers stud- 
ied this model or proposed other core-softened model po- 
tentials. 

Young and Alder carried out simulation of collapsing 
spheres (CS) system which is defined as U(r) = oo if 
r < a, U(r) = e if a < r < <j\ and U(r) =0 otherwise 
d, [H[ . Because of the computer power restrictions they 
managed to estimate the melting line of the system only 
in two-dimensional case. According to their's results the 
melting line demonstrates a maximum. They pointed 
out the similarity of the phase diagrams of the collapsing 
spheres system and experimental phase diagram of Ce 
and Cs. 

However, Young and Alder did not consider solid-solid 
phase transformation of collapsing spheres. This was 
done by Bolhuis and Frenkel who studied the CS system 
by Monte- Carlo simulation [f|. They considered isostruc- 
tural FCC-FCC transition in the collapsing spheres and 
square well systems. As it was shown in the cited paper, 
the isostructural transition takes place in both systems 
for sufficiently small soft core or square well widths. If 
the width of the soft core (well) increases the solid-solid 
transition disappears. 

Velasco et al. studied the collapsing spheres system 
in the frame of thermodynamic perturbation theory Q. 



Their's results are in good agrement with the Bolhuis and 
Frenkel publication [6j. They considered four values of 
the step width: <ti = 1.03; 1.08; 1.1 and 1.16. According 
to their's findings the first three systems demonstrate 
the isostructural solid-solid transformation with critical 
point while the last system does not have this transition. 
At the same time an FCC-BCC-FCC sequence occurs in 
the system with the largest step size. 

The qualitative phase behavior of collapsing spheres 
was studied by Stishov [8] . This article makes a qualita- 
tive estimation of the collapsing spheres phase diagram 
basing on the high- and low- temperature limits. It was 
shown that depending on the cti/ct ratio the phase dia- 
gram can change dramatically. This qualitative consider- 
ation is in good agreement with Monte-Carlo simulation 
Q and perturbation theory results [?J mentioned above. 
In Refs. 0, [Io| it was discussed the possibility of the 
liquid-liquid phase transition in CS system. 

One can use the geometrical reasons to identify the 
possible structures of the core softened systems. This 
method was applied by Jagla to a set of systems [ll[ • In 
this paper a parameter dependent potential was consid- 
ered. Depending on the value of the tuning parameter it 
can be continuously changed from a hard core plus linear 
ramp shape to hard core plus repulsive shoulder one. The 
minimal energy configurations corresponding to different 
pressures for the different potentials were monitored. It 
was established that depending on the pressure the num- 
ber of overlaps of the soft core changes producing a large 
number of possible geometrical configurations (Fig. 3 in 
11]). Basing on this argument one can expect a rich va- 



riety of structures in the core softened systems both in 
crystalline and disordered regions. 
Although the main focus of Ref. 



11] 



was on two- 
dimensional systems, the results were also generalized 
for three-dimensional case (Fig. 6 in [HI). These re- 
sults clearly state that the crystalline region of the linear 
ramp system demonstrates many phases with different 
symmetry. This fact identifies that core-softening is a 
very strong mechanism affecting the phase behavior of 
the system. 

The minimum energy configurations in the two- 
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dimensional collapsing spheres system was revised later 
in Refs. 



12 



13]. In this paper the MC simulation of the 
system was also carried out and complex phase behavior 
in liquid phase was observed. Depending on the density, 
the system can form different cluster phases. In particu- 
lar, long stripe clusters were observed at some densities 
and temperatures. 

The same system of two-dimensional collapsing spheres 
was studied further in Ref. [l4j ] . The stripe phases for 
a set of the step widths were observed. The authors 
pointed out that the stripe phase takes place only for 
some particular values of the step width. Although the 
paper was mainly focused on the two-dimensional system 
the authors claim that the conclusions are applicable for 
the three dimensional case too. 

Another soft core softened model in two dimensions 
was considered by Camp [15|, |l6j. Phase diagram ob- 



tained from NPT-MC was reported. This phase diagram 
consists from low- and high density solid phases, normal 
liquid and cluster liquid. The cluster liquid represents 
the stripe phase. 

Interestingly, the stripe phases were observed experi- 
mentally in the system of colloids in magnetic field 17 1. 



The configurations observed in this system closely resem- 
ble the ones obtained by Camp [ijj . This justifies that 
the core-softened models can be successfully applied to 
the real experimental systems. 

A strong method of the ground state phase calcula- 



tion was developed in Refs. [18j, [19]. This method allows 



to find the minimal energy configurations by applying 
generic algorithms carrying out the minimization. The 
collapsing sphere system with different step widths was 
considered in Refs. [3, [l{J. Three values of the step 
were considered: o\ — 1.5, <j\ = 4.5 and <j\ = 10.0. It 
was shown that the complexity of the system increases 
dramatically with increasing the step size. For example, 
for the system with a\ = 1.5 the number of phases re- 
ported in [3 is 6, for <j\ = 4.5 - 33 and for o\ = 10.0 - 
47. 

In our recent paper a smoothed repulsive shoulder sys- 
tem (RSS) was introduced [13]. This model is defined by 
the interparticle potential 
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(r ) W + \ £<yl ~ tanh ( fc °[ r ~ ai W> 



where fco = 10.0 and the step width a\ was varied. The 
dependence of the phase diagram on a\ was studied. The 
following values of a\ were considered: 1.15, 1.35, 1.55 
(Ref. [20]) and 1.80 (Ref. [ill]). In agreement with pre- 



vious discussion, the behavior of the system was found 
to be very complex. 

If the step width is relatively small then the crystalline 
region consists of low- and high density FCC phases sep- 
arated by a BCC region. However, already at o\ = 1.35 
the complexity of the phase diagram greatly increases. 
It is important to note that several crystalline structures 



appear in the system including not close packed struc- 
tures which can not be observed in LJ-like systems. Com- 
paring the phase diagrams for ai = 1.35, 1.55 and 1.8, 
one can see that the complexity of the phase diagram 
even increases with increasing the step width. 

Later the ground state of the RSS with a\ — 1.35 was 
studied in the work [22jj . The set of the phases reported 
in this work is different from the one in [20j. It can be 
related to the different choice of the phases the authors 
of these papers considered to be potentially stable. 

One can see that the behavior of the purely repulsive 
core-softened models is rather complex. However, the 
potentials of interaction between real particles have both 
repulsive and attractive parts. To take it into account 
one should add an attractive well to the repulsive core- 
softened model. The potentials of this type were reported 
as models of water, liquid metals, colloidal particles or 
polymer colloid mixtures and so on. This makes them 
important for the soft matter science and attracts much 
attention to such systems. 

The simplest core-softened system with attraction is 
the collapsing spheres system with additional attractive 
well (CSAW): 



co, r < (To 

£i, <Jq < r < (7i 

— £2, 0\ < r < (72 

0, r > cr 2 



(1) 



Liquid-liquid phase transition (LLPT) was found in 
such systems. LLPT in CSAW systems was intensively 
studied by molecular dynamics, integral equations theory 
and perturbation theory 2414291] . It was found that LLPT 
takes place only at some specific values of parameters 
£2/^1, ci/co an d 02 /oo- However, as far as we know, 
all papers on this topic suggest that LLPT is metastable 
with respect to crystallization. This makes extremely 
important to know also the melting line of the system. 
However, just a few publications consider the crystalline 
region of the phase diagram 3p|. 31 1. 

In the paper \^ the LLPT in CSAW was studied. 
The authors also investigated possible crystal structures 
in the system. They came to the conclusion that only 
one crystal phase is possible in the range of densities they 
considered. Taking into account huge variety of phases in 
the purely repulsive systems this conclusion looks rather 
surprising and makes it necessary to carry out further 
investigations of the influence of repulsive and attractive 
forces on the phase diagram of a model system. 

A detailed study of the phase diagrams of a family of 
core-softened systems with attractive well was carried out 
by Quigley and Probert 3(| 31 1. The potential used in 
their's works consisted of Lennard- Jones part plus a gaus- 
sian well. Choosing different parameters they changed 
the depth of the well and monitored the changes in the 
phase diagram. It was shown that FCC is not the only 
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number 
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(72 


Ao 


Ai 


A 2 


well depth 


1 


1.15 





0.5 


0.50 








2 


1.15 


1.35 


0.2 


0.5 


0.3 


0.4 


3 


1.15 


1.35 


0.07 


0.5 


0.43 


0.60 


4 


1.35 





0.5 


0.5 








5 


1.35 


1.80 


0.5 


0.60 


0.10 


0.20 


6 


1.35 


1.80 


0.5 


0.66 


0.16 


0.30 


7 


1.35 


1.80 


0.5 


0.7 


0.20 


0.4 



TABLE I: The potential parameters used in simulations (Eq. 
(2))- 



stable crystal structure for this kind of potentials. Simple 
hexagonal structure was also reported jsc" 



31|- 



The goal of this paper is to provide a detailed study 
of the influence of attractive and repulsive forces on the 
phase diagram of a core-softened system. For doing this 
we carry out the calculations of the phase diagrams of a 
set of systems which have the same functional form of the 
interatomic potential but parameters of this potential are 
different. In the present article we focus on the melting 
line of the system and the transformations between dif- 
ferent crystal phases. We also consider the influence of 
attraction on the thermodynamic anomalies in the core 
softened system. 



II. SYSTEM AND METHODS 

In the present study a system of particles interacting 
via the potential with "hard" core, repulsive shoulder and 
attractive well is investigated. This potential represents 
a generalization of our previous RSS model (UlEI and 
we call it repulsive shoulder system with attractive well 
(RSS-AW) potential. 

The general form of the potential is written as 



U(r) = e(-) 14 +A -Ai tanh(fci{r-c7i})+A 2 tanh(A; 2 {r-(7,}). 
r 

(2) 

Two sets of systems are considered: the ones with the 
step width equal to o\ = 1.15 and the ones with the step 
a% = 1.35. The parameters fci = fc 2 = 10.0 are fixed 
while parameters ai, (72, Ao, Ai and A 2 are varied to get 
the different potential shape. Seven sets of parameter are 
considered. They are summarized in Table 1. Figs. Q] 
and [2] show the potentials for the step width o\ = 1.15 
and <7i = 1.35 respectively. The parameters are chosen 
in such a way that the depth of attractive well becomes 
larger (see Table 1 and Figs. [T] and [2]). Below we 
denote the systems with different parameters as system 
1, system 2 and so on in accordance with the Table 1. 

In our previous publications [13, [HJ we discussed the 
phase diagrams of several purely repulsive systems, i.e. 
the systems with zero well depth. The complexity of 
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FIG. 1: (Color online) Family of the potentials with o\ = 1.15 
and different attractive wells. The curves are numerated in 
accordance with Table 1. 
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FIG. 2: (Color online) Family of the potentials with cri = 1.35 
and different attractive wells. The curves are numerated in 
accordance with Table 1. 



these phase diagrams was shown. Here we extend this 
study to the systems with attractive well and study the 
evolution of the system behavior with increasing attrac- 
tion. 

In order to get a hint about the possible crystal phases 
in the systems we measure the energies of ideal crys- 
tal lattices which correspond to zero temperature free 
energies. After that we simulate the most preferable 



4 



structures at finite temperature and monitor the limits 
of theirs stability by looking at the radial distribution 
functions. This procedure allows to find the approxi- 
mate regions of the thermodynamic parameters for free 
energy computations which are necessary for the exact 
determination of the transition lines. 

We simulate the system in NVT ensemble using 
Monte-Carlo method. The number of particles in the 
liquid or gas state simulation was set to 500 or 1000 and 
for crystal phases it varied between 500 and 1000 de- 
pending on the structure. The system was equilibrated 
for 10 6 MC step and the data were collected during 10 5 
MC steps. 

In order to find the transition points we carry out 
the free energy calculations for different phases and con- 
struct a common tangent to them. In our previous ar- 
ticle where we considered the purely repulsive poten- 
tials we computed the free energy of the liquid by in- 
tegrating the equation of state along an isotherm [35j |: 

F(.p)-Fid(p) _ 1 rf> P(p')-p'k B T , , j fh f 
Nk^T - fe^T Jo. 7 1 a P ■ 111 tJie CaSe 01 P°" 

tentials which contain an attractive part the situation 
is more complicated because of the possible gas - liq- 
uid transition. In order to avoid the difficulties con- 
nected to this transition we carry out calculation of 
free energies at high temperature above the gas - liq- 
uid critical point and then calculate the free energie s by 
integrating the internal energies along an isochor [35| : 

f V =J^(r,jy,v)4). 

Free energies of different crystal phases were deter- 
mined by the method of coupling to the Einstein crystal 

To improve the statistics (and to check for internal 
consistency) the free energy of the solid was computed 
at many dozens of different state points and fitted to 
multinomial function. The fitting function we used is 
a p,qT p V q , where T and V = X/p are the temperature and 
specific volume and powers p and q are related through 
V + 1 < N . The value N we used for the most of calcu- 
lations is 5. 

In this paper we use the dimensionless quantities: f = 
r/cr, P = Pa 3 /e, V = V/Na 3 = l/p,f = k B T/e. Since 
we use only these reduced units we omit the tilde marks. 



phase diagrams for the sake of comparison. As one can 
see from these figures, the phase diagram consists from 
liquid and three crystalline regions: low density FCC (ld- 
FCC), BCC and high density FCC (hd-FCC). Note that 
the shoulder o\ — 1.15 gives relatively small perturbation 
to the soft spheres potential $o( r ) = E (f) n with n = 14. 
As it was established in (3|| soft spheres freeze into BCC 
crystal if n < 6 and to FCC solid if n > 6. It means that 
the shoulder makes the potential effectively more " soft" . 

Figs. [3£b) and E£b) show the phase diagrams for 
the system 2, i.e. the system with a\ = 1.15 and the 
well depth w — 0.4. As it follows from the figures, the 
phase diagram is very similar. It again contains liquid 
region and three crystalline phases: ld-FCC, BCC and 
hd-FCC. One can see that the crystalline region shifts to 
the higher temperatures comparing to the previous case. 
This is clearly connected to the attractive forces which 
make the crystal phases more stable. 

The phase diagram for the system 3 is shown in Figs. 
Etc) and Etc). Here the attractive well is even deeper. 
It leads to further increase of stability of the crystalline 
phases. The melting line shifts to even higher tempera- 
tures. However, the general view of the phase diagram 
remains the same. 

Fig. [5] demonstrates all three phase diagrams in the 
same plot for the sake of comparison. Only P — T phase 
diagrams are shown since p— T plot would be overfilled by 
the data points. We observe several trends in the phase 
diagram evolution with increasing of the well depth. First 
of all, in the systems with attraction the ld-FCC phase is 
stable already at zero pressure while in the purely repul- 
sive system positive pressure is necessary to stabilize the 
crystal. Secondly the stability of the crystalline regions 
with respect to melting grows with increasing the attrac- 
tive well: melting temperature at the same density (or 
pressure) becomes higher as the attractive well becomes 
deeper. 

The main conclusion of this section is that in the case 
of the potentials 1 — 3 the attractive well shifts the phase 
diagram, but it does not change it's general shape. 

a\ = 1.35 (potentials 4-7) 



III. PHASE DIAGRAMS 
o i = 1.15 (potentials 1-3) 

In this section we consider the phase diagrams of the 
systems with step width equal to o\ = 1.15 and different 
depthes of the attractive well (systems 1 — 3 from Table 
1). These phase diagrams are shown in Figs. EKa) - 
(c) and H] (a) - (c). In our previous publication [20j we 
reported the phase diagrams for the purely repulsive po- 
tential (zero well depth). Figs. EJa) and Ufa) show these 



Now we turn to the potentials with larger width of 
the repulsive step o\ = 1.35. The phase diagram for the 
purely repulsive potential with o\ = 1.35 (potential 4) 
was published in our previous paper (2(j. After that the 
zero temperature phase diagram of this system was re- 
considered in the work Comparison of our phase 
diagram (Ref. [2(]j ]) and the zero temperature diagram 
from Ref. shows some differences. Body-centered 
orthorombic (BCO), body-centered tetrahedral (BCT) 
and (3 — Sn structures were reported to be stable into 
the region where we find face-centered tetrahedral (FCT) 
and simple cubic (SC) structures both at zero and finite 
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FIG. 3: (Color online) Phase diagrams of the systems with FIG. 4: (Color online) Phase diagrams of the systems with 
cri = 1.15 (system 1 - 3) in p — T coordinates, a - system 1, a\ = 1.15 (system 1 - 3) in P — T coordinates, a - system 1, 
b- system 2, c - system 3. b- system 2, c - system 3. 
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FIG. 5: (Color online) Phase diagrams of the systems 1-3 
in P — T coordinates. Empty symbols - system 1, semifilled 
symbols - system 2 and full symbols - system 3. 



temperature. Accordingly to our results FCT and SC 
phases are more stable then BCO and BCT in this re- 
gion, ft — Sn was not considered in our work. At the 
same time basing on the Ref. 



22j we find that simple 



hexagonal (SH) structure is stable in a certain range 
of densities. We carry out the free energy calculations 
for SH crystal and add the SH transition lines into the 
phase diagram (Figs. H](a) - (b)). 

The shape of this phase diagram can be intuitively ex- 
plained as follows. Since the potential is purely repulsive 
the overlap of both hard and soft cores is energetically 
unfavorable. It leads to the occurrence of the close packed 
structure at low densities (ld-FCC). However, when the 
pressure (or equivalently density) increases the particles 
penetrate the soft core. But since this core is repulsive, 
this penetration is energetically unfavorable and there- 
fore the system forms a set of structures with low co- 
ordination number - FCT (4 nearest neighbors), SC (6 
nearest neighbors) and SH (8 nearest neighbors). Fi- 
nally, the high density FCC (hd-FCC) crystal becomes 
stable which corresponds to the soft sphere limit. 

Taking this into account, one can expect that the ld- 
FCC phase should become more stable if the potential 
has also an attractive part in addition to the repulsive 
soft core. One can also expect that ld-FCC phase occurs 
al lower pressures comparable to the purely repulsive po- 
tential which corresponds to the arrangement of particles 
in the attractive well. 

Figs. [3(a) - (b) show the phase diagram for the system 
5, which corresponds to the well depth w = 0.2. 

Comparing the Figs. [5] and [51 one can see that the ld- 
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FIG. 6: (Color online) Phase diagram of the system 4 (Table 
1) in (a) p — T and (b) P — T coordinates. 



FCC phase do become more stable. While for the purely 
repulsive system the maximum temperature of ld-FCC 
phase is T max ^ = 0.19, for the system with attraction 
(system 5) the maximum temperature is T max ^ = 0.27. 
At the same time the region of stability of ld-FCC phase 
extends and goes into the negative pressure region. 

The FCT phase becomes stable in a wider density 
range. As a result the disordered gap presented in the 
system 4 phase diagram (Figs. [6] (a) - (b)) is filled by 
the FCT crystal. 

Another change in the phase diagram when the at- 
traction is added is that the SC crystal disappears. The 
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FIG. 7: (Color online) Phase diagram of the system 5 (Table 
1) in (a) p — T and (b) P — T coordinates. 



FIG. 8: (Color online) Phase diagram of the system 6 (Table 
1) in (a) p — T and (b) P — T coordinates. The inset of Fig. 
a shows the 1-FCC and FCT phases in large scale. 



region of SC crystal is now filled by FCT phase. The 
region of the SH phase is almost the same as in the case 
of small attraction. 

In order to trace the influence of the attraction on the 
phase diagram we consider the systems with deeper wells 
- w = 0.3 (system 6) and w — 0.4 (system 7). The phase 
diagrams of the system 6 are presented in the Figs. 8(a) 
and 8(b). 

One can see from Figs. |8](a) - (b) that the FCT phase 
goes to the lower temperatures region. Also a gap be- 
tween FCT and SH crystals appears again. One observes 
the triple point ld-FCC - FCT - Liquid at T = 0.11, but 
at lower temperatures the transition ld-FCC - FCT disap- 



pears. One can guess that FCT is not the most thermo- 
dynamically stable phase for this system. It means that 
if one finds a more stable phase all these contradictions 
will be solved. It gives the reasons to suggest that the 
larger attraction makes the system to form more complex 
crystal phases. The identification of these phases can be 
a very difficult task and goes beyond the scope of the 
present publication. 

The ld-FCC phase almost does not change. The maxi- 
mum on the melting line stays the same as in the previous 
case (potential 6). 

The range of pressures where SH structure is stable 
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FIG. 9: (Color online) Phase diagram of the system 7 (Table 
1) in (a) p — T and (b) P — T coordinates. 



reach its maximum, the temperature range of stability 
of this structure slightly decreases - from T^ ax 6 = 0.38 
to T^ ax 7 = 0.32. One can expect that, if the attraction 
becomes even stronger, the SH phase goes into the hd- 
FCC region and becomes metastable with respect to it. 
A large gap between ld-FCC and SH phases occurs. Most 
probably this gap is filled by some unknown structure as 
it was discussed above for the system 6. 

It is worth to trace the evolution of the Id- and hd- 
FCC phases with increasing attraction. From Figs. [TU] 
(a) - (b) one can see that when the attraction is added 
the ld-FCC phase becomes more stable with respect to 
temperature. However, increasing of the attraction does 
not stabilize the ld-FCC phase any more. At the same 
time the attraction moves the region of ld-FCC stability 
to the negative pressures. One can imagine that with fur- 
ther increase of attractive part the melting line will have 
negative slope even at zero pressure which qualitatively 
corresponds to the case of water. 

Attraction has much stronger effect on the hd-FCC 
melting line (Fig. [TT] (a) - (b)). Small attraction (w — 
0.2, system 5) almost does not affect the hd-FCC melting 
line, while increasing of the attraction moves it to the 
lower pressures. It makes the gap between ld-FCC and 
hd-FCC narrower. One can expect that at some value of 
the well depth the gap disappears and the crystal region 
of the phase diagram consists from two FCC phases only. 

Fig. [12] shows the specific volume jump over the 
crystal phase volume at hd-FCC melting for the systems 
4 — 7. Interestingly, the volume jump increases consider- 
ably with increasing the attractive well and reaches the 
value up to « 33%. Experimental volume jumps are usu- 
ally of the order of 5 — 10%, although for some substances 
it can reach much higher values. 

The example of the volume jump clearly shows the 
importance of the parameter £2/ei, where e\ is the step 
height and - well depth. Even at the same length 
parameters - <j\/<j and 02/ ' <J the character densities in 
the system can be strongly different. 



IV. THERMODYNAMIC ANOMALIES 



almost does not change again. This structure shows very 
weak dependence on the attraction. 

Figs. [9] (a) - (b) show the phase diagram of the sys- 
tem with well depth w = 0.4 (system 7). One can see 
from the figures that FCT phase disappears for this po- 
tential. The pressure range of the ld-FCC and SH struc- 
tures stay almost unchanged. At the same time the hd- 
FCC phase moves to the lower pressures region which 
makes the region of SH phase stability more narrow. In- 
terestingly the SH melting line looses the maximum ob- 
served for the previous systems. One can guess that this 
maximum goes under the hd-FCC melting line and be- 
comes metastable. As a result, since SH phase does not 



Finally we would like to discuss anomalous behavior of 
the core softened systems. In our previous works 2(| 21, 
[23l . |42| we showed that purely repulsive system (system 
4) demonstrates diffusion and density anomalies. Here 
we extend this study to the RSS-AW system. 

The problem of anomalous behavior of isotropic core 
softened system with attraction has already been dis- 
cussed by several authors 

In the recent publications [33, E(J the characteristic 
length and energies of the interatomic potential were 
fixed while the softness of the potential varied. It means 
that the focus of these publications is slightly different 
comparing with the present one. The authors found 
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FIG. 10: (Color online) ld-FCC phases for the systems 4-7 FIG. 11: (Color online) hd-FCC phases for the systems 4-7 
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anomalous behavior in the system they studied in a wide 
range of the potential parameters. 



In the publication [4l| the depth of attractive well is 
varied with all other parameters fixed. In this sense this 
work is very close to the present one. The deepest well 
considered in this work is —1.0 which is deeper then in 
our study. The main conclusion of the authors is that 
the anomalies shrink with increasing the well depth and 
finally disappear. 

Figs. [13] (a)-(c) present the diffusion anomaly for the 
systems 4, 5 and 7. One can see from these figures that 
the density range of anomalous diffusion does not change 
on increasing the well depth. At the same time the max- 



imum temperature of anomalous behavior increases. It 
means that small attraction stabilizes anomalous diffu- 
sion in the system. 

In order to see the origin of the difference in the qual- 
itative behavior of diffusion anomaly region between the 
present work and the publication 4l| we should notice 
that in our case the increasing of the well depth at the 
same time makes the potential softer (Fig. [2]). This cor- 
responds to making the parameter A, which controls the 
slope of the potential in the works [S, smaller. Ac- 
cordingly to these publications the anomalous region is 
wider for the smaller values of this parameter. Seems 
that this effect is larger in our case then shrinking of the 
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FIG. 12: (Color online) Volume jump over the crystal state 
specific volume for hd-FCC phase of the system 4 — 7 



anomalous region with increasing attraction. 

More detailed discussion of the anomalous behavior 
of the RSS-AW systems will be given in a subsequent 
publication. 



IV. CONCLUSIONS 

The present article shows the evolution of the phase 
diagrams of two core-softened systems with addition of 
attractive forces. It is shown that the general shape of 
the phase diagram and the set of preferable structures is 
mostly defined by repulsive shoulder while attractive well 
just shifts the phase diagram in p — T plane. Anomalous 
diffusion behavior is also considered. It is shown that 
the appearance of the diffusion anomaly if related to the 
soft-core but the small attraction stabilizes it at higher 
temperatures. 
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